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ABSTRACT 

Circumnuclear water masers can respond in two ways to changes in the luminosity 
of an active galactic nucleus. First, an increase in the X-ray luminosity can lead to an 
increase in the maser emissivity; and second, an increase in the intrinsic bolometric 
luminosity may result in a temporary decrease in the difference between the gas and 
dust temperature and a consequent decrease in the maser output. Whilst the latter 
effect can occur over a period shorter than the thermal timescale, the former effect 
cannot. Quantitative estimates of the response of the water maser emissivity to changes 
in either the X-ray or bolometric luminosity are presented, together with estimates of 
the relevant timescales. Either mechanism could account for recent observations by 
Gallimore et al. which suggest that the water maser variability in two widely separated 
regions of the circumnuclear gas in NGC 1068 have been coordinated by a signal from 
the active nucleus. For either mechanism, a minimum H2 density ~ 10^ cm~^ is needed 
to explain the observed variability timescale. 

Subject headings: masers - galaxies: active - galaxies: individual: NGC 1068 - 
galaxies: nuclei - molecular processes - radio lines: galaxies 

1. Introduction 

Luminous maser emissions in the 616 — 623 22 GHz line of water have now been detected from 
more than 20 active galaxies, a list of which is given in Table 2 of the recent review by Moran, 
Greenhill & Herrnstein (2000). In almost every case where the size of the emission region has been 
determined interferometrically, the water maser emission has been found to originate primarily 
within a parsec- or subparsec-sized region that is coincident with the active nucleus of the galaxy, 
suggesting that active galactic nuclei are the energy sources that ultimately power the observed 
water maser emission. In several sources, the geometry is strongly suggestive of a circumnuclear 
accretion disk (e.g. Miyoshi et al. 1995). 

Theoretical models for the collisional excitation of water masers (e.g. de Jong 1973, 
Neufeld & Melnick 1991) place requirements on the physical conditions in the emitting region. 
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Luminous water maser emission is favored by molecular hydrogen densities, n(H2), in the range 

cm ^, water abundances, a;(H20) = n(H20)/n(H2) ^ 10 ^, and gas temperatures 
^ 400 K. A minimum temperature requirement is imposed both by the excitation requirement 
for the maser transition and by the fact that water abundances are significantly enhanced at 
temperatures > 400 K by the reaction O + H2 ^ OH + H followed by OH + H2 ^ H2O + H (e.g. 
Elitzur & de Jong 1978; Neufeld, Lepp &; Melnick 1995). An additional requirement is that the gas 
temperature must significantly exceed the temperature of the radiation field at mid- to far-infrared 
wavelengths: thus some energy source that directly heats the circumnuclcar gas in active galaxies 
is needed to explain the observed water maser emission. X-rays from the active nucleus (Neufeld, 
Maloney & Conger 1994, hereafter NMC94), shock waves within a circumnuclear accretion disk 
(Maoz & McKee 1998), and viscous heating (Desch, Wallin & Watson 1998, hereafter DWW98) 
have all been proposed as the energy source responsible for heating the masing gas. 

Detailed models for water masers in X-ray heated gas have been presented by NMC94, who 
considered the water maser emission as a function of depth into an X-irradiated slab. Those 
models showed that the physical and chemical conditions within X-ray heated gas are conducive to 

the excitation of water masers, and that - given reasonable assumptions about the nature of the 
circumnuclear emission region - X-ray heating could account for the typical luminosities of water 
masers observed in AGN. CoUison and Watson (1995, hereafter CW95) subsequently made the 
important point that for certain parameters the maser emission can be significantly enhanced by 
the absorption of non- masing far- infrared water lines by dust, a process not included by NMC94. 

Since AGN are known to exhibit significant and rapid variability (e.g. Ulrich, Maraschi, &; 
Urry 1997) at optical, ultraviolet, and particularly X-ray wavelengths (e.g. Mushotzki, Done & 
Pounds 1993), maser models that invoke X-ray heating naturally raise the question of how the 
water maser emission would respond to changes in the incident X-ray and bolometric fluxes. This 
question, to be addressed this Letter, is also strongly motivated by observations of cooordinated 
maser variability in NGC 1068 that have been reported recently by Gallimore et al. (2000). In 
§2 below, I present the results of a general parameter study for the equilibrium maser emission 
from a warm dusty medium. In §3, I consider how the maser emission responds to changes in the 
illuminating X-ray flux or in the dust temperature within the masing medium. In §4, the results 
are discussed with speciflc reference to the water masers in NGC 1068. 

2. Water maser emission from a warm dusty medium 

The collisional excitation of water masers within a warm dusty medium has been modeled by 
Wallin &: Watson (1997) in the limit where the medium is optically thick to far- infrared continuum 
radiation. As previously pointed out by CW95, this limit often applies to circumnuclear gas in 
AGN, given the continuum dust opacities typically present (Draine &; Laor 1993). The effect of 
continuum opacity is to place a lower limit on the escape probability for non-masing far-infrared 
transitions of water, allowing maser emission to occur even at great depths within a circumnuclear 
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disk or gas cloud and enhancing the coohng rate due to far- infrared water emissions (DWW98). 

In the hmit where the medium is optically thick to far-infrared continuum radiation, the 
water level populations are determined solely by the local conditions and are dependent on the gas 

temperature, Tgas, the dust temperature, T^ustj the H2 density, n(H2), and the water abundance 
,t(H20) = n(H20)/n(H2). Since water abundances 10~^ are predicted for regions of luminous 
water mascr emission, I have constructed a grid of models for fixed x(H20) ~ 10"'^ that span a 
broad region in the three dimensional parameter space defined by Tgas, ^dust, and n(H2). 

The water level populations arc given by the usual equations of statistical equilibrium 
(e.g. CW95), with the rates for all radiative processes diminished by an escape probability 
pfj = [2ai,/(l + 2ay)](ln[e + 7r-V2ar.i])i/2, where aij is the ratio of dust continuum opacity to line 
center opacity for the transition from i to j (Hollenbach &; McKee 1979, CW95). For consistency 
with the previous work of CW95, we assume a continuum absorption coefficient of 10~^^nH cm^ 
for wavelengths A < 50 /um and 10^^^nH(50/im/A)'^ cm'^ for A > 50 /im, where ?t-h is the density 
of H nuclei. For comparison, the probability of escape from a semi-infinite plane parallel slab is 
Pij = (47r^/^Tij[lnT(j]^/^)~^ in the limit where the line center optical depth Tij » 1 (Hollenbach & 
McKee 1979). Comparing these expressions, we find that dust absorption dominates the effective 
escape of far-infrared water emissions for dust optical depths > (Stt^/^)"^ = 0.07, corresponding 
to column densities of H nuclei, A'^h ^ lO^^cm^^. Following Ncufcld & Kaufman (1993, hereafter 
NK93), I solved the equations of statistical equilibrium for 179 rotational states of ortho- water, 
adopting the collisional rate coefficients of Green, Maluendes &; McLean (1993) for transitions 
among the lowest 45 states and a simple extrapolation for transitions involving states of higher 
energy. The assumed line width, defined as velocity shift needed to reduce the opacity from its 
peak value by a factor e, was lkms~^. 

Once the statistical equilibrium equations have been solved, the cooling rate due to water 
emissions and the rate of maser emission under conditions of saturation may be obtained 
straightforwardly (e.g. Neufeld & Melnick 1991; NK93). The results I obtained for the maser 
emissivity are in good agreement with those presented previously by Wallin &: Watson (1997). If 
the gas is in thermal equilbrium, then the required heating rate is equal to the cooling rate, which 
- for the conditions of relevance here - is dominated by water cooling and gas-grain collisional 
cooling (Neufeld et al. 1995). Using the gas-grain cooling rate adopted by DWW98 and the water 
cooling rate obtained above, I computed the equilibrium heating rate, H, and used it in place of 
Tgas as one of the three dependent variables. 

Figure la shows a contour plot of the gas temperature, Tgas, and the maser efficiency, 
^sat — ^^^sat/l^-H-ff), as function of Tfjust and H/nu- Here $sat is the rate at which maser photons 
are emitted per unit volume, the superscript indicating that the results apply to the case where 
the dust optical depth is large (Watson &; Wallin 1997) and the subscript sat indicating that the 
maser emission is assumed to be saturated. The heating rate, H, is defined following the notation 
of Maloney, Hollenbach &; Tielens (1996, hereafter MHT96) as the heating rate per hydrogen 
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nucleus (with units erg s ^). The maser efficiency egg^^. is thus the fraction of the input power that 
emerges as 22 GHz maser radiation. The results shown in Figure la apply to a fixed molecular 
hydrogen density n(H2) = ^nn = 10^ cm~^. 

As required by the second law of thermodynamics, Figure 1 shows that Tgas — Tkust in the 
limit H ^ 0. There is a minimum value of the heating rate (defined by the curve e^^^^ = 0) 
below which Tgas — ?dust is insufficient to support a population inversion . The vertical line 
at H/uh ^ 1-2 X 10~^^ergcm^ s~^ represents the maximum X-ray heating rate that can be 
sustained without the gas being dissociated (NMC94). Thus the region for luminous maser action 
within an X-ray heated medium lies to the right of the curve e^^t = 0, to the left of the line 
H/riH = 1-2 X 10~^^ergcm^ s~^, and above the curve Tgas = 400 K (the latter requirement being 
needed for consistency with the large water abundance of 10^^ that has been assumed). Figure lb 
presents analogous results to those in Figure la, but now eg^t and Tgas are shown for a fixed dust 
temperature (500 K) as a function of n(H2) and H/tlh- 



3. Response of masers to a change in their environment 

Results shown in Figure 1 apply to conditions of thermal equilibrium. The effects of changing 
the heating rate and dust temperature can now be considered. The timescale for the gas to 
re-establish thermal equilibrium after a change in the heating rate or dust temperature is given 
by Tthcrm = lkT/2H = 0.027 (T/IOOOK) (n(H2)/109 cm'^)"! ([iJ/nHl/lO-^^ergcm^ s^^)- ^r. 
Let us first consider the effects of changing the heating rate. On timescales greater than Tthcrm) 
an increase in the heating rate increases the difference between the gas and dust temperatures 
and increases the rate of maser emission. The response is then conveniently characterized by 
the logarithmic derivative dln^^^^/dhiH = 1 -|- {dlne'^^^/dlnH), which is plotted in Figure 2 
(dashed curve) for n(H2) = 10^ cm~^ and T^ust = 500 K. Except for small values of H close to 
the threshold for maser action, the derivative (?ln<I>gat/91n/f lies within a factor of 2 of unity, 
implying that changes in the heating rate are accompanied by comparable fractional changes in 
the maser emission rate. 

The results plotted in Figure 2 are premised on assumptions of constant density and water 
abundance. The density is only expected to change on a dynamical timescale, which is typically 
much longer than rtherm- The molecular fraction, 2n(H2)/nH, only changes on the X-ray ionization 
timescale, which is given (MHT96) by ~ 4 (n(H2)/10'^ cm"^)"^ ([iJ/nnl/lC^^ergcm^ s~^)~Vr 
and is a factor ~ 200 greater than the cooling timescale. Even if the molecular fraction is constant, 
however, the water abundance can change in response to changes in the temperature. This effect, 
the inclusion of which lies beyond the scope of the present study, is unlikely to be important except 
when the gas temperature lies close to the threshold for rapid water production (Tgas ~ 400 K) or 
when the X-ray ionization rate lies close to the threshold for molecular dissociation. 

If the heating rate varies, the gas temperature and maser emissivity $gat cannot 
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respond on timescales any shorter than rtherm- However, $^3^^ can vary on timescales 

shorter than rtherm if the dust temperature changes rapidly. The thermal inertia of the 

dust grains is negligible, so the timescale on which the dust temperature can change 

is controlled by the diffusion time for the far-infrared continuum radiation[|, given by 

Trad ~ tfir£/c ~ 0.005[A''j//10^'* cm~^]^ [n(H2)/10^ cm^^J^^yr, where I is the distance below the 

irradiated surface of the cloud and tfir ~ 10 (A''///10^^ cm~^) (assumed > 1) is the corresponding 

optical depth at far-infrared wavelengths. 

Under circumstances where the dust temperature varies on timescales shorter than rtherm ~ 
as it might in response to variations in the incident bolometric flux from the AGN, for example 
- the gas temperature remains constant but the maser emissivity varies as a result of changes in 
the difference between the gas and dust temperatures, Tdifj = Tgas — Tdust- The response is now 
characterized by the derivative (91n$gj^^./91nrdust)Tgas) the subscript Tgas indicating that the gas 
temperature is held constant. The dotted curve in Figure 2 shows the value of this derivative, 
again for the case where n(H2) = 10^ cm^'^ and T^-ast = 500 K. Two features of this curve are 
particularly noteworthy. First, the value of (91n$gat/'91iiTdust)Tgas is always negative because Tdifj 
decreases as Tdust increases. Second, (91n<I>g3^^/91nTdust)Tgas can become much larger than unity 
for small values of H. This behavior occurs because the maser emissivity is primarily controlled by 
Tdifj; a small fractional change in T^ust causes a large fractional change in Tjig when T^\q <^ Tkust- 



4. Application to NGC 1068 

The results presented in Figures 1 and 2 can now be discussed with reference to recent 
observations of the Seyfert galaxy NGC 1068. Gallimore et al. (2000) have recently reported that 
two maser features in NGC 1068, symmetrically placed in Doppler velocity about the systemic 
velocity of the host galaxy, flared simultaneously with a rise time of at most 22 days. Given the 
emission region structure elucidated by earlier VLB A observations (Greenhill &; Gwinn 1997), 
the velocities of the maser features involved suggest that they are symmetrically located at 
distance d ~ 0.7 pc on either side of the nucleus. Unless the simultaneity of the flares is entirely 
happenstance, they were presumably coordinated by a signal that propagated from the nucleus. 

The results presented in Figure 2 suggest two possible mechanisms for how the maser flares 
might be coordinated: (1) variations in the central X-ray luminosity, Lx, might simultaneously 
change the heating rates in the masing regions; or (2) variations in the intrinsic bolometric 
luminosity, Lboi, might simultaneously change the dust temperatur^. Unfortunately, the average 



^This is the shortest timescale on which the local infrared radiation field can vary; note that the infrared flux as 
viewed from afar may have a characteristic variability timescale that is considerably longer than that of the local 
infrared radiation as a result of the different light travel times from the near and far sides of the irradiated disk or 
torus. 

^For the case of saturated masers, an third mechanism is possible: in principle, the observed flares might be 
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X-ray and bolometric luminosities are rather uncertain for NGC 1068, since the central source 
is obscured even for hard X-rays; I adopted the estimates of Pier et al. (1994), who compiled 
observations obtained over a wide range of wavelengths and assumed that 1% of the total intrinsic 
luminosity is reflected towards us. These estimates of the luminosities imply a 1 - 100 keV 
X-ray flux ~ 7 x 10^ ergcm"^ s~^ at the distance of the flaring masers and a characteristic 
dust temperature ~ (Lboi/47rci^cj)^/'^ ^ 500 K. Expressions given by MHT96 then imply an 
X-ray heating rate H ~ 5.6 x lO"^'^ (A^h/10^^ cm"^)""'^ erg s~^ at column density A'h below the 
irradiated surface of the maser emission region. For column densities greater than A^h ~ 10^^ cm~^, 
Thomson scattering becomes important and increases the effective column density appearing in 
the above equation for H hy a factor tt = 6.65 x 10~^^ A^h cm^. 

In Figure 3, relevant constraints are plotted]^ as a function of n(H.2) and A^h- Molecules are 
present above and to the right of the line H/rin = 1.2 x 10~^'^erg cm^ s~^. Maser action can occur 
below and to the left of the curve e^^^ = 0. The observed variability can be driven by changes 
in the X-ray heating rate for parameters to the left of the line Ttherm = 0-05 yr; or alternatively 
by changes in the dust temperature for the region to the left of the line Tj-ad = 0.05 yr and to 
the right of the line Ttherm = 0-05 yr. In either case, the minimum required II2 density for a 
variability timescale of 0.05 yr is ~ 108 cm ^. In Figure 4, the derivatives dln^^^^/dlnLx and 
(Sln'I'gg^^/SlnLboOrgas shown for several values of the density. Since the dust temperature 
increases as L^/^, the latter derivative is given by |(51n<I>sat/t^lii^dust)Tgas- 

Figures 3 and 4 indicate that the simultaneous flaring of the two maser features in NGC 1068 
could plausibly be coordinated either through variations in the heating rate or through variations 
in the dust temperature. The former mechanism would be unique to X-ray heating whereas 
the latter might apply in principle even if shock heating or viscous heating were responsible 
for the maser emission. Observationally, the mechanisms might be discriminated by searching 
for correlations between the maser emission and X-ray and bolometric luminosities, although in 
the case of NGC 1068 the fact that we observe only reflected radiation from the nucleus will 
introduce uncertain delays that will be difficult to remove. Maser variability induced by a varying 
X-ray heating rate will be characterized by a positive correlation with Lx (with a time delay 
determined by the different light travel times, of course), while variability induced by a varying 



coordinated by a varying 22 GHz continuum flux from a compact central source. Since the maser radiation from 
the flaring regions travels tangentially toward us, the flux we receive from those regions might increase if the radial 
flux of seed radiation from the central source diminished. This explanation, however, would require that the maser 
emission is controlled by seed radiation from a compact continuum source; in NGC 1068 the observed distribution of 
maser emission close to the systemic velocity argues against that being the case. 

^The results plotted here apply to a gas cloud that is illuminated at normal incidence such that the column 
densities for X-ray shielding and for the escape of far-infrared radiation are identical. For oblique illumination at 
angle 9 to the normal, the column density A''h refers to the shielding column density, the timescale r, ad is reduced by 
a factor sec^B, and the dust temperature should be reduced by a factor (sec^)^''*. An additional assumption is that 
dust absorption dominates the removal of far-infrared line photons, requiring A^h cosO ^ 10^^ cm^^. 
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dust temperature will be characterized by a negative correlation with Lboi- 
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Fig. 1. — Maser efficiency (solid contours labeled with 10^ x e^^^) and gas temperature (dashed 
contours labeled with Tgas in K). Results are shown as a function of heating rate, H, and either 
dust temperature (Figure la, top panel, for which the H2 density is fixed at 10^ cm""^) or H2 density 
(Figure lb, bottom panel, for which the dust temperature is fixed at 500 K). The assumed water 
abundance for both panels is 10~^ relative to H2 . The region favorable to efficient maser emission 
is the unshaded region where the gas temperature is large enough ( ^ 400 K) to yield a substantial 
water abundance and where the heating rate supports a value of (Tgas — Tdust) large enough to 
permit maser action but is not so large as to dissociate molecules. 
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logjo H/n„ (erg cm^ s ') 

Fig. 2. — Response of maser emissivity to small changes in the heating rate and dust temperature, 
for a dust temperature of 500 K and a H2 density of 10^ cm^'^. Values of the derivatives 
dln^^^^/dlnH (dashed curve) and (51n$ga^^/91nTdust)Tgas (dotted curve) are shown as a function 
of the heating rate, H (see text). The solid curve shows the maser emissivity (value plotted is 

10^ X 4f ) 



- 11 - 




Fig. 3. — Parameter space for maser action and variability in clouds 0.7 pc distant from the active 
nucleus of NGC 1068. Results are shown as a function of the H2 density and the column density 

for X-ray shielding. Regions of parameter space for which there is no maser emission are shaded. 
The dashed hnc shows where the thermal timcscalc for the gas is 0.05 yr, and the dotted line shows 
where the radiative diffusion timescale governing changes in the dust temperature is 0.05 yr. The 
coordinated maser variability reported by Gallimore et al. (2000) can be driven by changes in the 
X-ray heating rate for parameters to the left of the line rtherm = 0.05 yr; or alternatively by changes 
in the dust temperature for the region to the left of the line Trad = 0.05 yr and to the right of the 
line rtherm = 0.05 yr. 



- 12 - 




logio N„ (cm-^) 

Fig. 4. — Same as Figure 2, but now plotted as a function of the column density for X-ray shielding. 
Results apply to a cloud of H2 density 10^ cm~^ that is 0.7 pc distant from the active nucleus of 
NGC 1068. Arrows indicate the range of column densities for which the thermal timescale for the 
gas and the radiative diffusion timescale governing changes in the dust temperature are smaller 
than 0.05 yr. 



